
Observations of the Crab pulsar with the MAGIC telescopes 
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Abstract: We report on the observations of the Crab pulsar with the MAGIC telesopes. Data were taken both in the 
mono-mode (> 25 GeV) and in the stereo-mode (> 50 GeV). Clear signals from the two peaks were detected with both 
modes and the phase resolved energy spectra were calculated. By comparing with the measurements done by Fermi-LAT, 
we found that the energy spectra of the Crab pulsar does not follow a power law with an exponential cutoff, but that it 
extends as a power law after the break at around 5 GeV. This suggests that the emission above 25 GeV is not dominated 
by the curvatura radiation, which is inconsistent with the standard prediction of the OG and SG models. 
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1 Introduction 

The Crab nebula is the compact object left over after a his- 
toric supernova explosion that occurred in the year 1054 
A.D.. The pulsar B0531+21 (also commonly named Crab 
pulsar) is located at its center, and emits strong pulsed ra- 
diation in a wide energy range from radio to high energy 
gamma-rays. 

The Crab pulsar and a few other pulsars are among the 
brightest known sources at 1 GeV. However, a spectral 
steepening made their detection above 10 GeV elusive de- 
spite numerous efforts. The energy thresholds of imag- 
ing atmospheric Cherenkov telescopes (IACTs) were, in 
general, too high, while the gamma-ray collection area of 
satellite-borne detectors were too small to detect pulsars 
above 10 GeV. On the other hand, a precise measurement 
of the energy spectrum at and above the steepening leads to 
an important verification for the standard pulsar models. In 
the case of the Polar Cap (PC) model Q , so-called super- 
exponential cutoff is expected, while the Outer Gap (OG) 
model |2] and Slot Gap (SG) model Q, predict a clear ex- 
ponential cutoff. The highest energy of the detected pho- 
tons can be directly converted to the lower limit on the dis- 
tance of the emission region from the stellar surface, which 
should be a few times the stellar radius according to the PC 
model. 

In 2008, the MAGIC telescope detected the Crab pulsar 
above 25 GeV |4] with the newly implemented trigger 
system, the Sum trigger (5). This detection excluded 
the PC model. In August 2008, the new satellite borne 



gamma-ray detector with 1 m 2 collection area, Fermi-LAT, 
became operational and it could measure the spectra of 
gamma-ray pulsars up to a few tens of GeV. The spectra 
measured by Fermi-LAT could be described with a power 
law with an exponential cutoff, which also rejected the 
polar cap model and supported the OG and the SG model. 

However, the cutoff energy of the Crab pulsar spectrum 
determined by Fermi-LAT was ~ 6 GeV, while MAGIC 
detected the signal above 25 GeV. In order to verify the 
exponential cutoff spectrum, i.e. OG and SG models, the 
precise comparison of the energy spectra measured by the 
two instruments is needed. Here we present the spectral 
study of Crab pulsar, using the public Fermi-LAT data and 
four years of MAGIC data recorded by the single telescope 
and the stereoscopic system. 

2 MAGIC telescope 

The MAGIC telescope is a new generation IACT located 
on the Canary island of La Palma (27.8° N, 17.8° W, 2225 
m asl). It consists of two telescopes with a reflector diam- 
eter of 17 m. The first telescope was build in 2002-2003 
and have been operational since 2004. Thanks to the world 
largest reflector, the energy threshold of the first MAGIC 
telescope with the standard trigger is 60 GeV, that is the 
lowest among IACTs. In order to detect gamma-ray pul- 
sars, the new trigger system called Sum trigger was devel- 
oped and implemented in October 2007. It reduced the en- 
ergy threshold further, down to 25 GeV, which resulted in 
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the detection of the Crab pulsar [4]. In 2009, the second 
MAGIC telescope was build ~ 80 m apart from the first 
telescope. The second one is basically a clone of the first 
one, while the Sum trigger system is not yet implemented 
to it. 

We observed the Crab pulsar with the stereoscopic mode 
from 2009. The stereo trigger requires a coincidence of the 
triggers of both telescopes. For a technical reason, the Sum 
trigger in the first telescope cannot participate in the stereo 
trigger, i.e., stereoscopic observations were based on the 
standard trigger for both telescopes. The energy threshold 
of the stereo mode is about 50 GeV. 

3 Mono-mode observations 

MAGIC observed the Crab pulsar with a single telescope 
with the Sum trigger in winter 2007-2008 and winter 2008- 
2009. After the careful data selection, total effective ob- 
servation time was 25 hours and 34 hours for the first and 
the second campaign, respectively. The energy threshold 
of these observations are 25 GeV. Normally, IACT tech- 
nique utilizes many image parameters to distinguish be- 
tween hadron events and gamma-ray events. However, in 
the case of mono-mode observations at the very low en- 
ergy regime below 60 GeV, the image parameter are almost 
powerless except for the Hillas parameter ALPHA. There- 
fore, the hadron background rejection was done only based 
on ALPHA. 

The light curve of the Crab pulsar obtained with the mono- 
mode observation is shown in the upper panel of Fig. [T] 
Following the usual convention [6] of PI (phase interval 
-0.06 to 0.04) and P2 (0.32 to 0.43), the numbers of ex- 
cess events in PI and P2 are 6200 ± 1400 (4.3 a) and 
11300±1500(7.4a). By summing up PI and P2, the excess 
corresponds to 7.5a. The background level was estimated 
by using the so-called Off Pulse phase (0.52 - 0.88) 0. 

Based on these excess events, the phase resolved energy 
spectra of the Crab pulsar above 25 GeV were computed 
as shown in Fig. [2] They can be well described by power 
laws and the best fit parameters are summarized in Table 
Q] The energy spectrum measured by Fermi-LAT is also 
shown in the same figure. For the Fermi-LAT points, 1 
year of the Fermi-LAT data (from August 2008 to August 
2009) were used. The best fit parameters obtained by the 
unbinned likelihood analysis are summarized in Table [T] 
The continuation from the Fermi-LAT measurements to the 
MAGIC measurements is rather smooth, while it is clear 
that the exponential cutoff spectrum determined by Fermi- 
LAT is not valid in MAGIC energies. A detail statistical 
analysis showed the inconsistency amounts to 6.7a, 3.0a, 
and 5.8a for PI + P2, PI and P2, respectively. Further 
details of the results of the mono-mode observations are 
presented in 0, @ and (TD|. 



4 Stereo-mode observations 

As mentioned before, MAGIC observed the Crab pul- 
sar with the stereo-mode since 2009. Though the energy 
threshold of the stereo-mode (50 GeV) is higher than the 
mono-observation with the Sum trigger (25 GeV), the VHE 
tail of the spectrum allows us to detect the Crab pulsar 
above 50 GeV. The advantage of the stereo-mode obser- 
vation is higher background rejection power, better angu- 
lar resolution and better energy resolution. In the case of 
the stereo-mode observations, by using images from both 
telescopes, one can reconstruct the arrival direction of the 
gamma-rays better than in an ALPHA analysis. The esti- 
mation of the shower maximum height is also more precise 
than in the mono-mode, which leads to the higher back- 
ground rejection and the better energy resolution. 

The lower panel of Fig. Q] shows the light curve of the Crab 
pulsar obtained by stereo-mode observations. The total ob- 
servation time used in this analysis is 73 hours. Since the 
pulses are much narrower than the conventionally defined 
PI and P2 phases, the evaluation of the statistical signif- 
icance of the excess was done by Z^Q-test, H-test and \ 2 
test, which gave 8.6a, 6.4a and 7.7a, respectively. By fit- 
ting two Gaussians to the two peaks, the peak positions 
are estimated to be 0.005 ± 0.003 and 0.3996 ± 0.0014, 
while the corresponding FWHMs are 0.025 ± 0.007 and 
0.026 ± 0.004. By defining the signal phases as ±2a of 
the fitted Gaussian around the peaks, the significance of 
the excess is 10.4a, 5.5a and 9.9 a for PI + P2, PI and P2, 
respectively. The phase resolved energy spectra are also 
calculated and shown in Fig. [3] The dark red squares de- 
note the spectra of the ±2a phase interval described above, 
while the yellow squares denote the ones with the conven- 
tional P1/P2 definitions. They smoothly connect with the 
mono-mode measurements and follow a power law. Fur- 
ther details of the results of the stereo-mode observations 
are presented in ifTTl . 

5 Discussion and Conclusion 

The new measurements with MAGIC-mono and MAGIC- 
stereo discovered that the energy spectrum of the Crab pul- 
sar does not roll off as fast as the exponential cutoff but 
extends as a power law after the break at around 5 GeV. 
This suggests that the emission above 25 GeV is not dom- 
inated by the curvature radiation as the standard OG and 
SG models predict. Further theoretical studies and more 
observations are needed to understand the emission mech- 
anism of the Crab pulsar. A theoretical interpretation of the 
MAGIC measurements by K. Hirotani is presented in J9), 
|[T0l and ED. 
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Figure 1 : The light curve of the Crab pulsar obtained by 
MAGIC above 25 GeV with the mono-mode (top) and 
above 50 GeV with the stereo mode (bottom). The effec- 
tive observation time is 59 hours for mono and 73 hours for 
stereo. 
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Figure 2: The phase resolved energy spectra of the Crab pulsar measured by Fermi-LAT and MAGIC with the Sum trigger 
(mono-mode). The energy spectra measured by Fermi-LAT is consistent with a power law with an exponential cutoff, 
while the MAGIC measurements above 25 GeV deviate from their extrapolation. The inconsistencies amount to 6.7a, 
3. Oct and 5. 8 a level for PI + P2, PI and P2, respectively. 





Fermi-LAT 


MAGIC-Mono 


Phase 


h [10- iU cm- 2 s- i MeV- i ] 


Ti 


E c 


/so lO-^nrt^TeV- 1 ] 


r 2 


PI +P2 


1.94 ±0.05 


1.98 ±0.02 


4.5 ±0.3 


14.9 ± 2.9 


3.4 ± 0.5 


PI 


1.29 ±0.04 


1.99 ±0.02 


3.7 ±0.3 


4.5 ±2.3 


3.1 ± 1.0 


P2 


0.67 ±0.02 


1.95 ±0.03 


5.9 ±0.7 


10.0 ±1.9 


3.4 ± 0.5 



Table 1: The best fit parameters of the spectra for different phase intervals. The second to forth columns are obtained 
with Fermi-LAT data assuming the spectral shape of = h(E/l [GeV])- Fl exp(-E/E c ), while the fifth and sixth 

columns are obtained with MAGIC-mono data assuming the spectral shape of = f 3 o(E/30 [GeV])- r2 




Figure 3: The phase resolved energy spectrum of the Crab pulsar measured by Fermi-LAT, MAGIC-mono and MAGIC- 
stereo. After the break at around 5 GeV, the spectra follow power laws. (For the difference between the yellow points and 
dark red points, see text.) 



